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OO We present a family of solvable models of interacting particles in high di- 

o 
1^ 

o 

Qn^ glassy regime with aging effects. The interaction is controlled by a parameter 



I 

o 



mensionalities without quenched disorder. We show that the models have a 



C^ p. For p = 2 we obtain matrix models and for p > 2 'tensor' models. We 

concentrate on the cases p = 2 which we study analytically and numerically. 
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Disordered systems such as structurally disordered materials (glasses) and spin-glasses 
are characterized by complex free-energy landscapes and non-equilibrium phenomena with 
time-scales exceeding the experimental times [1-5]. The glass transition being mainly kinetic, 
one has to face the full dynamical problem in order to grasp many aspects of the physical 
behavior of these systems. This is a very difficult problem for realistic models. 

As regards to spin-glass systems, it has been recently noticed that simple yet microscop- 
ically inspired models, for which the mean-field solution is exact [4], have both solvable and 
non-trivial dynamics which agrees at least qualitatively with that of realistic systems [6-8]. 
Spin-glasses are described by Hamiltonians with quenched (random) disorder included by 
hand and do not have an ordered ground state to which the system can possibly decay. 
Structural glasses instead have an ordered ground-state though when cooled sufficiently fast 
they remain in a regime of (self-induced) disorder. Various aspects of the glassy behavior 
have been studied with phenomenological scenarios, dynamical theories, toy kinetic models 
and numerical simulations (see e.g. refs. [1-3]). 

However, unlike the case of spin-glasses, there has been a lack of simple, solvable Hamil- 
tonians capturing the essentials of structural glasses. Several recent papers investigate mean- 
field models which despite having no explicit quenched disorder still have a glassy transition 
[9,10,12]. In particular, the model of Ref. [12] has been shown to obey dynamical equations 
identical to those of the p-spin glass solved in Ref. [6]. It is not clear whether the emergence 
of such glassy behaviour in these models is a consequence of using 'complicated' [i.e. quasi 
random) Hamiltonians. The object of this note is to present a family of very simple models 
that can be interpreted as interacting particles in high dimensionalities and are suitable to 
describe the supercooled liquid phase and the glassy regime. 

We consider a system of n particles moving in an iV- dimensional space, in the limit 
n,iV — )• oo [13]. In terms of the coordinates s" of the a-th particle the Hamiltonian reads: 

H^-Y.{s^-s'r, (1) 

where s" = (s",...,s^) and a,h = l,...,n. We consider the case (it) in which the particles 



are constrained to move on a spherical surface | s"" |^= N Va, and the case in which they 
can occupy only the vertices of a hypercube s^ = ±1. The constraints counterbalance 
the repulsive force and confine the particles. The range of the interaction is smaller the 
larger the value of p. We choose n such that the force term f3'V°'E on each particle a has 
components of order one and involves all other particles if they are distributed randomly 



(i.e. if (5" • s^) = (9(v iV)). This implies the set of choices f3 ~ 0{J ^'^_i ) and, in particular, 
we use n = 2aNP-'^/p\ with jS ~ 0(1) or n = 2aN/p\ with jS ~ 0{N'^-p^^), both of which 
yield non-trivial models. 

In this paper we concentrate in the model with p = 2: 



Z = e^ Tr{s} exp 



-M{S^SS^S) 



(2) 



aN 

where S is the N X n rectangular matrix of elements 5", and Tr^^} runs over either the 
spherical or the Ising measure. The cases a > 1 and a < 1 are related; if a > 1 the free- 
energy satisfies F[a,f3) = 1/a F[l/a,f3 = af3) + (a — l)iV^. We also report very briefly 
some results for the model with p = 3,4 and, furthermore, for the it version of (2) but with 
symmetrical S. 

The partition function (2) is that of a rectangular matrix model. In the spherical case it 
can be solved exactly; it does not have a glassy phase transition but behaves like a 'liquid' 
at all temperatures. The it model has, instead, three phases: 

i. A high-temperature ('liquid') phase in which its free-energy coincides (up to a con- 
stant) with that of the spherical model. 

ii. Low-temperature T < Tf ('crystalline') states that dominate the low-temperature 
Gibbs-measure. The statical transition is first order. 

Hi. A dynamical low temperature ('glass') regime T < Tg. 

Our strategy is as follows: we first show the equivalence between spherical and it models 
in the high-temperature phase, and we solve the spherical model for all temperatures. Then 
we briefly discuss the nature of the 'crystalline' states of the it model. Finally, we verify 



numerically, using Montecarlo dynamics, these analytical results. From the simulations it 
becomes clear that the it model has a dynamical transition at a temperature Tg-. below Tg 
a large enough system asymptotically reaches an energy that is higher than the analytical 
continuation of the high-temperature phase, and it never (in the limit N — )• oo) 'falls' into a 
'crystalline' state. Furthermore, we show that the dynamics in the 'glass' phase is remarkably 
similar to that of both mean-field and realistic spin-glasses. We expect that similar results 
should be obtained with e.g. Langevin or Glauber dynamics. 

Let us first perform a Gaussian decoupling of the quadratic term in Eq. (2) by introducing 
auxiliary variables Mab- 

■ IP 

a,b 

The measure DM depends on the M ensemble we consider, both choices DM 



Z = x{f3) I DM Tr{,.} exp I "^V ^ E ^abs" ■ sA . (3) 



UaKbdMab JiV/47rexp(-iV/4 trM^) [M real, symmetric with zero diagonal) and DM 



Ua<bdRe{Mab) <jN/ATrUa<bdIm{Mab) <J N / Att exp{- N / A tiM^) {M hermitian) are possi- 
ble. Then % = 1 and % = e^^ , respectively. 
For the spherical model we have 



's.H = x(/3) / n ^ exp [^ E z"l / DM exp 



-^trln \dSabz'' + 2iJ^Mab 
I \ a \ a , 



(4) 



Since z" are only iV variables and the exponent is symmetrical w.r.t. the permutations of the 
indices of z", an admissible saddle point is —iz°' \s.p.= Z^, Va. This implies that the solution 
coincides, in the large iV limit, with that of the same model with the global constraint: 
EW^r = nN. 

The model written in this form, with hermitian M, can be recognized as a form of the 
generalized Penner model [14] or of the Penner-Kontsevich model before integrating over z" 
[15]. 



We now prove that the high-temperature expansions for the free-energies of the spherical 
and the it models coincide to all orders in (3. Let us go back to the it version of Eq. (3), 
and add a term = | Y^a ^a{N— \ s"' p). We later fix Ua appropiately. With the symmetric 
decoupling Eq. (3) becomes: 



Z± = e^^c'"^ [dm exp 



N In Tr^^a} exp ( - 2 XI ^"'D^bs'' 

\ ab 



(5) 



with Dab = "^iy f3 / aMab + Sabi^a- Performing the trace over the 5" after a Gaussian decoupling 
the exponent in (5), NW{M,u), is 

W{M,u) = -hv\n{D + I) + n\n2 + \nlY[Xaexp(-^J2^-^~'-bXb-J2^+ - ) 

O' \ ab ct / 

(6) 

where A~^ = D~^ + I and the dots stand for the higher orders of the series of ln(cos Aq). 

The off-diagonal elements of the propagator Aa^ are (9(l/v iV). We now choose the Ua 
to impose that the diagonal elements are also of that order. This yields for Ua the implicit 
relations [D + I]~^ = 1, and W{M,u{M)) becomes 

N 1 

W{M, z/) = - ^ z/„ - -tr In (D + 7) + n In 2 + 'diag'(M) (7) 

The term 'diag', stands for the diagrams generated by the higher orders in A^ in (6) with the 
propagator A. These are functions of M which can then be treated as perturbations of tr M^ 
in the integral over M (for example, the first term is 'Yabi-^ab]'^)- They are in turn vertices 
of the 'fat diagrams' (with propagator tr M"^) which lead to non-planar diagrams, and hence 
can be neglected in the large N limit (order by order in f3). This is in contradistinction to 
the vertices in the 'fat diagrams' coming from the trace of the logarithm in (7) (of the form 
tr M'^) which cannot be neglected for any k. 

Identifying Ua <-> if3za/a — 1, the condition for the u"" now coincides with the saddle point 
equation for the z" in the spherical model. Hence both free energies coincide to leading order 
in iV, up to an additive constant in the entropy. 



In order to solve the spherical model we return to Eq. (2). Assuming a > 1 and using 
the large N equivalence with the global constraint we get, up to constant multiplicative 
factors: 

Zsph ~ e^""" I'^t' jli ^^^ ^^P (-^'^ - /^^[^^] ) • (8) 

The Xiy N are the 'diagonal' values of S in its canonical form and E[xi\ reads 

We define a symmetric density of eigenvalues p(^x) = | Y^i[S{x — Xi) + S[x + Xi)] to get the 
equation (in the continuous limit) 



{a-l)= fdy-^^, (10) 

J X — V 



^-^-(a-l)=/cZy^M 

g g J X -y 



where g = 2aT. Following Refs. [16,17] we propose the generating function 



F{z) = If! _ /^ _ (^J) _ If J{z^-a^){z^-P){z' + c') . (11) 

g g z g 

The jump across the cut of F yields p. The constants a, 6, c and the saddle point value 
of fi are determined by the conditions: i. p ^ 1/z as z — )• oo (normalization), ii. F[z) 
has no pole in z = and Hi. the constraint X^ a;? = n = aN. Defining a± = a^ it 6^, 
these conditions yield fifth degree polynomials for a_|_, a^, c^ and fi. The energy density 
is obtained from a contour integral of z^Fi^z) around the cut and gives e = E/[nN) = 
(l/a)[-l + al/{32ga){al + 4a^ + Sa+c^)]. In fig. 1 we plot e vs. T for a = 1. 

At T ~ 0, E = N^[{a - 1) + T/A] for a > 1 (and E = N^T/{Aa^) for a < 1), showing 
explicitly that the entropy tends to — oo logarithmically at T = 0. Then the entropy of the 
high-temperature solution vanishes at some temperature To > (for a = 1, To ~ 0.15) and 
there must be a temperature Tf > To in which the it system undergoes a static transition. 
Two possibilities then arise: either the low-T partition function is dominated by a glassy 
state or by an ordered 'crystalline state'. For the spherical model a little algebra shows that 
the ground states are obtained by constructing the N X n matrix S with rows (columns) of 



orthogonal vectors of norm yN {\/n) for n > N [N > n). For the it model, we are not 
always able to construct such a matrix with ±1 elements, though for some values of n,iV 
this is possible {e.g. n = N = 2'^), in which case the energy of the lowest configuration 
coincides with that of the spherical model. For these iV, n these states dominate the low-T 
partition function; for other values we have observed numerically other states which though 
non-optimal play the same role at low-T. 

Since these states have very small entropy, Tf is approximately the temperature at which 
the free-energy of the high-temperature coincides with their energy (for a = 1, T/ ~ 0.41). 

The important point is that these 'crystalline' ground states are different from the glassy 
regime and are not reached dynamically by large enough systems, as we shall see. 

In fig. 1 we show the analytic and numerical results for the temperature dependence of 
the energy. The numerical points were obtained by means of a slow annealing in temperature, 
the energy densities obtained in this way ('dynamical energies') tend to have a well-defined 
limit with iV, for N large enough. For the spherical model the 'dynamical energy' and the 
analytically predicted static energy coincide for all temperatures. In the high temperature 
phase the 'dynamical energies' of the spherical and ±1 models also coincide. 

Below the glass temperature Tg ~ .5 (> T/), the 'dynamical energy' of the it model is 
higher than that of the spherical model. For the cases in which the ground state is known 
{e.g. N = n = 2^) we have checked that the system takes a time that grows fast with N to 
'fall' into a crystalline state. 

In order to undertand the nature of the 'glass' phase, we have studied the decay of the 
auto-correlation functions C{t,t') = ^ I {Nn)(^a s°'{t) ■ s°'{t')). Above Tg for the it model, 
and at all temperatures for the spherical model, the averaged auto-correlations decay with 
a time translationally-invariant law {C{t + t^^t^) = C{t)) as in a system in equilibrium. 
Below Tg the it model presents aging effects [5]: the auto-correlation functions shown in 
fig. 2 are clearly those of a system away from equilibrium even for very long times. The 



auto-correlations C{t + t^^t^) have, for large t^t^, an approximate tjl^ dependence (for 
large tjl^ C ~ [t/tw)~'°^) and are strikingly similar to the corresponding ones of mean-field 
[6-8] and d = 3 spin-glasses [18]. 

Let us mention in passing that we have performed dynamical simulations of the it model 
(1) with p = 3 (n ~ iV^, /3 ~ 0(1)), with p = 4 (n ~ iV, /3 ~ 0{1/N)) and also of (2) with 
S^ = S. In all these cases we have found glassy behavior. 

In order to make contact with the theory of systems with quenched disorder we first note 
that the if in the system described by the Eq. (1) we study the dynamics of one particle, 
while fixing the other n — 1 particles in random positions, we obtain neural network models 
(at negative temperature) [19]; systems known to have a spin-glass behaviour. 

A different strategy, that can be used to solve all the models we have described, consists 
in proposing a model with quenched disorder which is expected to have the same dynamical 
behavior [9,10]. For the case p = 2 the natural choice is to consider the Hamiltonian (1) 
with the s" defined by s" = U{ai,bj) ^^, where t^ are either spherically constrained or ±1. 
The tensor U{ai,bj) is a quenched random variable taken from the ensemble of orthogonal 
transformations "R"-^^ -^ -j^nxN ^^_^_ preserving the norm || t ||^= I](^")^). One then 
considers (quenched) averaged quantities over U. 

The analysis of the statics uses the replica trick and goes along the lines of Ref. [10]. We 
do not describe it in detail here, but report some results: If the t^ are spherically constrained 
the disordered model coincides with the non-disordered spherical model for all temperatures. 
If the t'^ take values ±1, the high-temperature phase (uncoupled replicas) also coincides with 
the spherical model. In order to analize the low-temperature (glassy) phase of the disordered 
model one ends up with a replicated matrix model. 

Perhaps more ambitiously, one could try to face the full dynamical problem as developed 
for mean-field spin glasses [6]. 
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Fig. 1 Numerical values for the energy density vs. temperature for the spherical (crosses) 
and lb (diamonds) models. N = n = 50, a = 1. The solid line shows the analytical solution 
for the spherical model. 
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Fig. 2 Aging curves for the ± model for a = 1 at T = 0.2 (< Tg); In — In plot, C{t + 1^, t^) 
vs. (t/t^) for t^ = 1000,300,100,30, t up to 3000. N = n = 500. 
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